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We study the non-linear dynamics of quantum fields in matter and radiation 
dominated universes, using the non-equilibrium field theory approach com- 



<D , bined with the non-perturbative Hartree and the large N approximations. 



We examine the phenomenon of explosive particle production due to spin- 



odal instabilities and parametric amplification in expanding universes with 
■ and without symmetry breaking. For a variety of initial conditions, we com- 

pute the evolution of the inflaton, its quantum fluctuations, and the equation 
of state. We find explosive growth of quantum fluctuations, although par- 
ticle production is somewhat sensitive to the expansion of the universe. In 
the large N limit for symmetry breaking scenarios, we determine generic late 
time solutions for any flat Friedman-Robertson- Walker cosmology. We also 
present a complete and numerically implementable renormalization scheme 
for the equation of motion and the energy momentum tensor in flat FRW 
cosmologies. In this scheme the renormalization constants are independent of 
time and of the initial conditions. 
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I. INTRODUCTION 



Over the last fifteen years, there has been a sustained effort to study the evolution of 
scalar fields within the context of the early universe. These efforts have been largely fueled 
by the introduction of inflation which has been shown to be a possible solution of the 
horizon and flatness problems. In addition to model building and other discussions of how 
scalar field theories can produce inflation, the question of the reheating of the universe, the 
transfer of energy from the inflaton to other particle modes, has received much attention 
||. The reason for this is that exponential expansion during inflation causes the universe to 
become very cold, as the energy in fields other than the inflaton is redshifted away. In order 
to achieve the standard results of nucleosynthesis, and possibly other early processes, the 
universe must be reheated to temperatures above those at which these important processes 
take place. Early efforts to account for the necessary reheating introduced ad hoc decay 
widths to the inflaton, assuming the energy transfer occurred through single particle decay 
Q]. However, more recently, it has been realized that there are much more efficient processes, 
those of either spinodal decomposition in the case of a (new) inflationary phase transition 
and parametric amplification in the case of chaotic inflationary scenarios [^H7|,|To Hl6| , p!9 



Such mechanisms, in which the rate of energy transfer grows exponentially, are referred to 
generally as preheating. 

Recently we have treated the dynamical processes of inflation and preheating [|7| within 
a fully non-equilibrium formalism [p0 -|2~4}l. Such a treatment is necessary due to the ex- 
pansion of the universe and the nature of the rapidly evolving dynamics of preheating. 
Non-equilibrium analyses have shown a wealth of new phenomena which were missed in 
equilibrium studies [[^,p|,p!6|-[T9|| . 

It has also become clear that the early analyses in which the inflaton is treated as a 
classical field or within perturbation theory are, in many cases, inadequate. In fact, despite 
the tiny couplings usually assumed for the inflaton, the unstable growth of modes during 
preheating causes the dynamics to become non-perturbative f25H27|. 

In this article, we study the dynamics of scalar fields in an expanding isotropic universe 
using the non-equilibrium closed time path formalism (CTP) pOfl , keeping track of the evo- 
lution of both the zero mode of the inflaton and its fluctuations. We treat the dynamics 
using two non-perturbative schemes. These are the Hartree approximation appropriate to 
a theory with discrete symmetry and the leading order large N approximation of an O(N) 
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vector model which describes theories with continuous symmetry, satisfying the correspond- 
ing Ward identities. . Both of these are mean field theory approximations, and 
as such they cannot account for the particle scattering processes that would allow the uni- 
verse to reenter the hot Big Bang scenario after inflation. Eventually, these processes should 
be taken into account to determine the final reheating temperature. Here we will concern 
ourselves only with processes occurring before thermalization. 

In particular, we will study the process of preheating using a wide range of initial condi- 
tions while the study of the inflationary stage will be discussed elsewhere [0. Taking our 
cues from the evolution of the equation of state in the inflationary and preheating phases, 
we analyze the dynamics of preheating in fixed radiation and matter dominated Friedman- 
Robert son- Walker cosmologies both analytically and numerically. We follow the equation of 
state during the evolution to ensure that our evolution obeys the appropriate gravitational 
dynamics. 

In the next section, we set up the closed time path formalism, describe our models and 
the two approximations we will use to study the evolution, and write down our evolution 
equations for the zero mode and the fluctuations. In section III, we discuss some important 
issues regarding the renormalization aspects of the problem, including the renormalization 
of the energy and pressure densities in an expanding background. Section IV begins with 
early time solutions in the slow roll scenario followed by a full numerical analysis of the 
various cases. We conclude the section with late time analyses of the large N evolution 
equations in the case of a symmetry broken potential. In the conclusions, we contrast our 
work with other analyses in the literature and discuss avenues which may be pursued to 
further improve our knowledge of these important cosmological problems. An appendix is 
provided with a discussion of our choice of initial conditions and their physical implications. 

Our main results are as follows. In the situations we analyze, we find that the expansion of 
the universe allows for significant particle production, although this production is somewhat 
sensitive to the exact expansion rate and is effectively shut off for high enough rates. In the 
case of a symmetry broken potential, we determine that in the large N limit, the quantum 
fluctuations decrease for late times as l/a 2 (t), while these fluctuations and the zero mode 
satisfy a sum rule consistent with Goldstone's theorem. In addition to these results, we 
present a consistent renormalization of the energy momentum tensor in a flat FRW spacetime 
within both the large N and the Hartree approximations. Such a result is an essential 
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component of any consistent analysis of the backreaction problem in an expanding universe. 

We compute the renormalized energy density, e, and the pressure, p, as a function of 
time. Averaging over the field oscillations, we find immediately after preheating a cold 
matter equation of state (p = 0) in the slow roll scenarios. In chaotic scenarios, the equation 
of state just after preheating is between that of radiation (p = e/3) and matter where 
the matter dominated regime is reached only for late times. The timescale over which the 
equation of state becomes matter dominated depends on the distribution of created particles 
in momentum space in addition to the approximation scheme implemented. 



II. THE FORMALISM AND MODELS 



We present here the framework of the non-equilibrium closed time path formalism. For 
a more complete discussion, the reader is referred to 0, or the alternative approach given 

0- 



in 



The time evolution of a system is determined in the Schrodinger picture by the functional 
Liouville equation 

i^- = [H(t), P (t)\, (2.1) 

where p is the density matrix and we allow for an explicitly time dependent Hamiltonian as is 
necessary to treat quantum fields in a time dependent background. Formally, the solutions to 
this equation for the time evolving density matrix are given by the time evolution operator, 
U(t,t'), in the form 

p(t) = U(t,t )p(t )U- 1 (t 1 t ). (2.2) 

The quantity p(to) determines the initial condition for the evolution. We choose this initial 
condition to describe a state of local equilibrium in conformal time, which is also identified 
with the conformal adiabatic vacuum for short wavelengths. In the appendix we provide an 
analysis and discussion of different initial conditions and their physical content within the 
context of expanding cosmologies. 



Given the evolution of the density matrix (|2.2|), ensemble averages of operators are given 



by the expression (again in the Schrodinger picture) 



m)) = Tr[U(to t t)OU(t,t')U(t',to)p(t )] 

Trp(t ) 
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where we have inserted the identity, U(t,t')U(t ,t) with t' an arbitrary time which will be 
taken to infinity. The state is first evolved forward from the initial time to to t when the 
operator is inserted. We then evolve this state forward to time t' and back again to the 
initial time. 

The actual evolution of various quantities in the theory can now be evaluated by either 
constructing the appropriate Green functions as in , or by choosing an explicit ansatz for 
the functional form of the time dependent density matrix so that the trace in ( |2.3| ) may 
be explicitly evaluated as a functional integral (see [fL7|l ). The methods are equivalent, and 
provide the results which will be presented below for the cases of interest. 

Since we are currently interested in the problem of preheating at the end of inflation, we 
will work in a spatially flat Friedman-Robertson- Walker background with scale factor a(t) 
and line element: 

ds 2 = dt 2 - a 2 (t) dx 2 . (2.4) 
Our Lagrangian density has the form 



V M $V$ - V($) 



(2.5) 



A. Hartree Approximation 

In the Hartree approximation, our theory is that of a single component scalar field, 
$(x, t), with the Z 2 symmetry $ — > — $. The potential can be written as: 

V($) = i(m 2 + ^)$ 2 + ^$ 4 , (2.6) 

where 1Z is the Ricci scalar. We decompose the field into its zero mode, (pit) = ($(x, £)), 
and fluctuations i/i(x,t) about it: 

$(f,t) = <f>{t) +ip(x,t). (2.7) 

The potential ( |2.6| ) may then be expanded in terms of these fields. 

The Hartree approximation is achieved by making the potential quadratic in the fluctu- 
ation field ip by invoking the factorization 

7/> 3 (f,t)->3(V 2 (x,t))V(f,t), (2.8) 
i) 4 (x, t) -> 6(^ 2 (f, t))ip 2 (x, t) - 3(^ 2 (x, t)) 2 . (2.9) 
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This factorization yields a quadratic theory in which the effects of interactions are encoded 
in the time dependent mass which is determined self-consistently. 

The equations of motion for the zero mode and the fluctuations are given by the tadpole 
equation 



Introducing the Fourier mode functions, Uk(t), they can be written as: 

k*) + 3^|| Mt) + (m 2 + (R(t)) m + j <t,\t) + j = 



^ +3 ls + 4 +m2+W)+ ^ 2(t) n w ' 2(i)> 



2(2vr) ; 



(2.10) 

(2.11) 
(2.12) 
(2.13) 



The initial conditions on the mode functions are 



U k (t ) 



1 



with the frequencies u) k {to) given by 



Uk(to 



a(t c 



a(t c 



U k (t 



Wk(to)= k 2 + M 2 (t ) 



1/2 



M 2 (t) = a 2 (t) 



(2.14) 



m 2 + {i-l/Q)n(t) + ^ 2 {t) + \{V{t)) 



(2.15) 



A detailed analysis and discussion of the choice of initial conditions and the frequencies 
( P-15| ) is provided in the appendix. As discussed there, this choice corresponds to the 
large-fc modes being in the conformal adiabatic vacuum state. In what follows we will 
subtract the composite operator ip 2 {t) at the initial time and absorb the term |(V ;2 (to)) in 
a renormalization of the mass. Furthermore, we choose the scale of time such that a(t ) = 1 
in both radiation and matter dominated cosmologies. 

B. Large N Limit 



To discuss the large N limit, we now treat $ as an iV-component vector in a theory with 
the continuous O(N) symmetry. The potential is 



6 



V($) = \{m 2 + ■ $ + A($ . (2.16) 

We now break up the field $ into a single scalar field a and an iV — 1 component vector 
7? as $ = (a, ff) and allow the a field to have a non-zero expectation value. Taking the 
decomposition 



<r(2,t) = VN<p(t) + x&t), (2.17) 

where y/N(j)(t) is the expectation value of a and x is the fluctuation about this value. If we 
write the "pion" field as 

AT-l 



7t(x,t) = ^(x,t)(l,l,...,l), (2.18) 
we reach the leading order large N limit by assuming the factorization 

^ 4 (f,t) - 2(^ 2 (f,t))^ 2 (f,t) - (^ 4 (£,t)>, (2.19) 

with 

(^ 2 > = 0(1), (x 2 ) = 0(l), = 0(1). (2.20) 

We see that since there are N —1 pion fields, contributions from the field x can be neglected 
in the formal limit as they are of order 1/N with respect those of ip and <fi. 



Again, we determine the equations of motion via the condition (|2.10| ). The zero mode 
equation becomes 

fa) + 3 ^§<^) + t™ 2 + + + \ m (^ 2 (t)) = o, (2.2i) 

while the equations for the modes, the fluctuation, and the initial conditions are identical 
to the Hartree case given by equations (|2.12| ) - (|2.14 ). Notice that we have used identical 



notations in the two cases to avoid cluttering and also to stress the similarity between the two 
approximations. In particular, we note that the only difference in the expressions for the two 
cases |eqs.( [2.12"D and ( |2.21|) , respectively] is a factor of three appearing in the self interaction 



term in the equations for the zero mode. However, as we will see, the two approximations 
are describing theories with distinct symmetries and there will be qualitative differences in 
the results. 

An important point to note in the large iV equations of motion is that the form of the 
equation for the zero mode (|2.21| ) is the same as for the k = mode function ( |2.12| ). It will 



be this identity that allows solutions of these equations in a symmetry broken scenario to 
satisfy Goldstone's theorem. 



III. RENORMALIZATION AND ENERGY MOMENTUM 



Upon examination of the equal time correlator ( 2.13|) , one finds that the integral is 
divergent and thus must be regulated. This can be done by a number of methods, but 
we require a scheme which is amenable to numerical calculation. We therefore introduce a 
large momentum cutoff which renders the integral finite, and one finds that it is possible 
to remove the terms depending both quadratically and logarithmically on the cutoff by a 
renormalization of the parameters of the theory [|TT | . After these terms are subtracted, the 
cutoff may be taken to infinity, and the remaining quantity is both physical and finite. A 
similar process is required to regulate the expressions for the energy density and the pressure 



as will be described in more detail below [18 



In terms of the variables introduced in section II above, the renormalization of (|2.13|) 
proceeds almost identically in the Hartree and large N approximations. The WKB analysis 
that reveals the large-/c behavior of the mode functions is described in detail in the appendix 
wherein we quote the relevant expressions for the large-/c behavior for the mode functions 
and their derivatives. Denoting bare and renormalized parameters with subscripts b and r 
respectively, and the momentum cutoff and subtraction point as A and k respectively, and 
using the results of the appendix, we find the following renormalization scheme: 



m 2 b + 



X b A 2 
I671- 2 a 2 (t) 



1 + 



A;, 



167T 



■ln(A/«) 



A; 



1 

£-4 



7 A r ln(A/V)/167r 2: 



16tt [ 



•(£ P -l/6)ln(A/«) 



(3.1) 
(3.2) 
(3.3) 



where 7 is a combinatorial factor taking on the value 7 = 1 in the large N limit and 7 = 1/3 
in the Hartree approximation. The subtracted equal time correlator is now given by 

^ d 3 k i\U k (t)\ 2 1 



+ 



(2tt) 3 
B{k-K) 
4k 3 



2ka?(t) 



A 



l/6)K(t) + Tn* + f(<j?(t) + (rp(t)) r 



(3.4) 



A further finite subtraction at to is performed and absorbed in a further finite and 
time independent renormalization of the mass. Notice in particular that in order for the 
renormalization of the mass to be time independent, we must require that the cutoff A be 
fixed in physical coordinates and therefore have the form A cx a(t). 
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Our treatment of the renormalization of the energy momentum tensor is similar to the 



approach of |30|], extended to the non-perturbative Hartree and large iV approximations. 
The expressions for the expectation values of the energy density, e, and the trace of the 
stress energy, e — 3p, where p is the pressure density are: 



e 
N 



1 :. 



H — m 
2 



2/2 



+ 



7A 



+ 



m 



27A a 



e — 3p 
N 



a 

7 A 



(3.5) 



-</> z + 2m 2 (b 2 4 



+ 



2m 4 
7A 



+ r + 3- 



60(^ 2 ) + 



((W) 2 ) + (2 - 60m 2 (^ 2 ) - OT(1 - 60<tf 2 > 



A , 



+ -[(2 - Qorm + (rr - etmmr), 



(3.6) 



where we have used the equations of motion in deriving this expression for the trace ( |3.6|) . 
The quantities = —TZ and G® = — 3(d/a) 2 are the trace and the time-time components 



of the Einstein curvature tensor, (if) 2 ) is given by equation ( 2.13|) , (ip 2 ) r by 
have defined the following integrals: 

d 3 k 



<(W) 2 > 
<^ 2 > 



2(2vr) 3 

d 3 k 
2(2tt) 3 



k 2 \U k {t)Y 



u k {t)f 



and we 

(3.7) 
(3.8) 



The composite operator (ipip) is symmetrized by removing a normal ordering constant to 
yield 



d 3 k d\U k (t)\' 



dt 



(3.9) 



Each of these integrals is divergent and must be regularized. We proceed in the same manner 
as above, imposing an ultraviolet cutoff, A, and computing the high k expansions of the £4, 
this time to fourth order in 1/k. We find the following divergences in e and e — 3p: 



A 4 



+ 



A 2 



div 



167r 2 a 4 167r 2 a 2 



2(e-l/6)G° + m 2 + ^(0 2 + (^ 2 ) ? 



+ 



ln(A/«) 

167T 2 



- - 2 y (0 2 + <V> 2 >r) - |(0 2 + <^ 2 >,) 2 
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+ 2(£ r - l/6)Gg (m 2 + ^(0 2 + ^ 2 > r )^) + (£ r - l/6) 2 tf ° 



3p\ 



AT 



div 



+ 



6(£ r 

A 2 
167r 2 a 2 
ln(A//c) 

16vr 2 



l/6)-^( 
7 7 a 2 dt 



2 

^ 2 + <n 



(3.10) 



2(^-1/6)^+12(^-1/6)- 



2m 2 + A r (0 2 + (^ 2 ) r ) 



-2m* - 2m 2 A r ( 



A 



+ 2(£ r - 1/6)G£ (m 2 + 



)) + & - 1/6) 2 #; 



- 6(e r - 1/6) 



A r d 2 
Yrft 2 ' 



+ (^ 2 )r)+3 



a X r d 



(3.11) 



a 2 dt 

The quantities ifg and H£ are the time-time component and trace of a geometrical tensor 
given by the variation with respect to the metric of higher derivative terms appearing in the 
gravitational action (such as 1Z 2 ). For the present case, they are given in terms of the Ricci 
scalar by the expressions 



a 2 



a* 
d- 

3-n 

a 



12 



(3.12) 
(3.13) 



Eqs.( |3.10"|) and ( |3.11| ) are closely related to eqs.(3.17) in ref. [fl8l| . However, they are not 
identical since the initial conditions chosen in ref. were different from the ones selected 
in the present paper. Furthermore, we had £ r = in ref. [jn| . 

The energy momentum is made finite by subtraction of the divergent pieces ( |3.10| ) and 
( [3. 11| ) from the expressions for the energy density (|3.5| ) and the trace ( |3.6| ). Within the 
context of covariant regularization schemes such as dimensional regularization and covariant 
point splitting [29|, such a procedure has been shown to adequately renormalize couplings 
appearing in the semi-classical Einstein's equation, extended to account for higher derivative 
terms and a possible cosmological constant, K. This equation has the form 



+ am + 



Kgl 



(3.14) 



where Gat is Newton's constant and a is the coupling to the higher order gravitational 
term. However, regularization via an ultraviolet cutoff is not a covariant scheme and we 
find that the quadratic and quartic divergence structure of the energy momentum in this 
scheme does not have the correct form to consistently renormalize the parameters of the 
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semi-classical theory. Nevertheless, the subtraction procedure described here produces a 
stress energy tensor which is both finite and covariantly conserved in addition to being 
amenable to numerical study. An alternative renormalized computational scheme in which 



covariant regularization is possible is presented by Baacke, Heitmann, and Patzold |31|| . It 
should be possible to extend such a scheme to expanding spacetimes. 

IV. EVOLUTION 

We focus our study of the evolution on radiation or matter dominated cosmologies, as 
the case for de Sitter expansion has been studied previously fl8 j. We write the scalefactor as 



a(t) = (t/to) n with n = 1/2 and n = 2/3 corresponding to radiation and matter dominated 
backgrounds respectively. Note that the value of to determines the initial Hubble constant 
since 

rr ( , \ a(*o) n 

H ( t o) = -7TT = r- 
a(t ) t 

We now solve the system of equations ( ^.11| ) - ( |2.13| ) in the Hartree approximation, with 
( f2.21| ) replacing ( |2.11| ) in the large N limit. We begin by presenting an early time analysis 



of the slow roll scenario. We then undertake a thorough numerical investigation of various 
cases of interest. For the symmetry broken case, we also provide an investigation of the late 
time behavior of the zero mode and the quantum fluctuations. 

In what follows, we scale all variables in terms of the magnitude of the renormalized 
mass, taking \mf \ = 1. We also define the variable 

v 2 (t) = x<f> 2 (t)/2 

and write 

gZ(t) = X(^(t)) r /2 , g = A/Stt 2 = 1(T 12 . 

We drop the subscript r denoting the renormalized parameters, and we will assume minimal 
coupling to the curvature, £ r = 0. In the most of the cases of interest, 71 <C 1, so that finite 
£ r has little effect. 

A. Early Time Solutions for Slow Roll 



For early times in a slow roll scenario [m 2 = —1, rj(ta) <ti 1], we can neglect in eqs.( p.lil) 



or (|2.21|) and in eq.( 2.12 ) both the quadratic and cubic terms in T](t) as well as the quantum 
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fluctuations (ip 2 (t)) r [recall that (?/> 2 (to))r = 0]. Thus, the differential equations for the zero 
mode ( 2. 11| ) or (|2.21|) and the mode functions (|2.12|) become linear equations. In terms of 
the scaled variables introduced above, with a(t) = (t/t ) n (n = 2/3 for a matter dominated 
cosmology while n — 1/2 for a radiation dominated cosmology) we have: 



3t^ 

v(t) + T v(t) - v(t) 



(4.1) 



k 2 



d 2 3n d 

dl? + TJt + {t/t ) 2n 



- i 



U k (t) = . 



The solutions to the zero mode equation ( [4.1[ ) are 

T 1 {t) = cr v i v {t) + dr v K u (t) 



(4.2) 



(4.3) 



where v = (3n — l)/2, and I v {t) and K u {t) are modified Bessel functions. The coefficients, c 
and d, are determined by the initial conditions on r\. For r](t ) = i] and fj(t ) = 0, we have: 



c = Vo C 1 



V 



K v (t Q ) - -K v (t ) 
to 



d 



I v (t ) - -I, (t ) 



(4.4) 
(4.5) 



Taking the asymptotic forms of the modified Bessel functions, we find that for intermediate 
times rj(t) grows as 



t»i c 



^7T 



t -3n/2 e t 



9n 2 — Qn _ / 1 



(4.6) 



We see that r/(t) grows very quickly in time, and the approximations ( |4.1| ) and ( |4.2| ) will 
quickly break down. For the case shown in figure 1 (with n = 2/3, H(to) = 0.1, T](to) = 10~ 7 , 
and 7)(to) = 0), we find that this approximation is valid up to t — to — 10. 

The equations for the mode functions ( |4.2|) can be solved in closed form for the modes 
in the case of a radiation dominated cosmology with n = 1/2. The solutions are 



U k (t) = c k e-tu(\-^\ 



Z,2t ] +d k e-*M\ ^~ 



3 k 2 t% n 3 

— , — , 2t 
I 2 '2' 



(4.7) 



Here, U(-) and M(-) are confluent hypergeometric functions |32[ (in another common nota- 
tion, M(-) = 1 F 1 (-)), and the c k and d k are coefficients determined by the initial conditions 
fl2.14j ) on the modes. The solutions can also be written in terms of parabolic cylinder func- 
tions. 
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For large t we have the asymptotic form 
U k (t) ^ d k e *(2t)-( 3 / 4 +W72) 



np(3 (hi™))- 

Z 1 \4 2 

+ c k e-* (2t)(- 3 / 4 +( fc *o) 2 /2) ^ l + P^j 



1 + 0(1 



(4.8) 



Again, these expressions only apply for intermediate times before the nonlinearities have 
grown significantly 

B. Numerical Analysis 

We now present the numerical analysis of the dynamical evolution of scalar fields in time 
dependent, matter and radiation dominated cosmo logical backgrounds. We use initial values 
of the Hubble constant such that H(t Q ) > 0.1. For expansion rates much less than this value 
the evolution will look similar to Minkowski space, which has been studied in great detail 
elsewhere |7]|y33]. As will be seen, the equation of state found numerically is, in the majority 



of cases, that of cold matter. We therefore use matter dominated expansion for the evolution 
in much of the analysis that follows. While it presents some inconsistency at late times, the 
evolution in radiation dominated universes remains largely unchanged, although there is 
greater initial growth of quantum fluctuations due to the scale factor growing more slowly 
in time. Using the large N and Hartree approximations to study theories with continuous 
and discrete symmetries respectively, we treat three important cases. They are 1) m? < 0, 
T](t ) < 1; 2) m 2 < 0, T](t ) > 1; 3) m 2 > 0, rj(t ) > 1. 

In presenting the figures, we have shifted the origin of time such that t — > t' = t — to- 
This places the initial time, to, at the origin. In these shifted coordinates, the scalefactor is 
given by 

a(t) = 



T 

where, once again, n = 2/3 and n = 1/2 in matter and radiation dominated backgrounds 
respectively, and the value of r is determined by the Hubble constant at the initial time: 

ff(to = 0) = ~. 

T 

Case 1: m 2 < 0, i](to) 1. This is the case of an early universe phase transition in 
which there is little or no biasing in the initial configuration (by biasing we mean that the 
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initial conditions break the 77 — > —rj symmetry). The transition occurs from an initial tem- 
perature above the critical temperature, T > T c , which is quenched at t to the temperature 
Tf <C T c . This change in temperature due to the rapid expansion of the universe is modeled 
here by an instantaneous change in the mass from an initial value m 2 = T 2 /T 2 — 1 to a final 
value m 2 = —1. We will use the value m 2 = 1 in what follows. This quench approximation 
is necessary since the low momentum frequencies (|2.15|) appearing in our initial conditions 
fl2.14j ) are complex for negative mass squared and small t]{to). An alternative choice is to 
use initial frequencies given by 



k z + M ^(to)tanh 



'k 2 + M 2 (t )^ 1/2 



\M 2 (t )\ t 

These frequencies have the attractive feature that they match the conformal adiabatic fre- 
quencies given by ( |2.15| ) for large values of k while remaining positive for small k. We find 
that such a choice of initial conditions changes the quantitative value of the particle number 
by a few percent, but leaves the qualitative results unchanged. 

While this case should show some qualitative features of the corresponding process of 
preheating in slow-roll inflation, we note that since quantum fluctuations can grow to be 
large during the de Sitter phase (see |]18[), a proper treatment of preheating after slow- roll 
inflation must account for the full gravitational backreaction. This will be the subject of a 
future article |34|. In the present case, we impose a background cosmology dominated by 



ordinary radiation or matter. We plot the the zero mode t]{t), the equal time correlator 
g£(t), the total number of produced particles gN(t) (see the appendix for a discussion of 
our definition of particles), the number of particles gNk(t) as a function of wavenumber for 
both intermediate and late times, and the ratio of the pressure and energy densities p(t)/e(t) 
(giving the equation of state). 

Figure la-e shows these quantities in the large N approximation for a matter dominated 
cosmology with an initial condition on the zero mode given by 77(^0 = 0) = 10~ 7 , r)(t = 0) = 
and for an initial expansion rate of H(to) = 0.1. This choice for the initial value of 77 stems 
from the fact that the quantum fluctuations only have time to grow significantly for initial 
values satisfying rj(t ) <C ^fg\ for values r](t ) ^> ^fg the evolution is essentially classical. 
This result is clear from the intermediate time dependence of the zero mode and the low 
momentum mode functions given by the expressions ([4.6| ) and ( |4.9|) respectively. 

After the initial growth of the fluctuation gE (fig. lb) we see that the zero mode (fig. 
la) approaches the value given by the minimum of the tree level potential, 77 = 1, while #E 
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decays for late times as 

9S c 



a 2 it) t 4 / 3 ' 

For these late times, the Ward identity corresponding to the O(N) symmetry of the field 
theory is satisfied, enforcing the condition 

-l + V 2 (t)+gi:(t) = 0. (4.9) 

Hence, the zero mode approaches the classical minimum as 

" 2(i) - 1 - im ■ 

Figure lc depicts the number of particles produced. After an initial burst of particle 
production, the number of particles settles down to a relatively constant value. Notice that 
the number of particles produced is approximately of order 1/g. In figure Id, we show the 
number of particles as a function of the wavenumber, k. For intermediate times we see the 
simple structure depicted by the dashed line in the figure, while for late times this quantity 
becomes concentrated more at low values of the momentum k. 

Finally, figure le shows that the field begins with a de Sitter equation of state p = —e but 
evolves quickly to a state dominated by ordinary matter, with an equation of state (averaged 
over the oscillation timescale) p = 0. This last result is a bit surprising as one expects from 
the condition ( |4.9|) that the particles produced in the final state are massless Goldstone 
bosons (pions) which should have the equation of state of radiation. However, as shown in 
figure Id, the produced particles are of low momentum, k <^1, and while the effective mass 
of the particles is zero to very high accuracy when averaged over the oscillation timescale, 
the effective mass makes small oscillations about zero so that the dispersion relation for 
these particles differs from that of radiation. In addition, since the produced particles have 
little energy, the contribution to the energy density from the zero mode, which contributes 
to a cold matter equation of state, remains significant. 

In figures 2a-e we show the same situation depicted in figure 1 using the Hartree ap- 
proximation. The initial condition on the zero mode is 77 (to = 0) = \f?> ■ 10~ 7 ; the factor 
of y/3 appears due to the different scaling in the zero mode equations, (|2.11|) and fl2.21Q , 



which causes the minimum of the tree level effective potential in the Hartree approximation 
to have a value of 77 = a/3. Again, the Hubble constant has the value H(t ) = 0.1. Here, 
we see again that there is an initial burst of particle production as gE (fig. 2b) grows large. 
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However, the zero mode (fig. 2a) quickly reaches the minimum of the potential and the 
condition 

- 1 + 7] 2 (t)/3 + gZ(t) = (4.10) 

is approximately satisfied by forcing the value of gE quickly to zero. There are somewhat 
fewer particles produced here compared to the large N case, and the distribution of particles 
is more extended. Since the effective mass of the particles is nonzero, we expect a matter 
dominated equation of state (fig 2e) for later times. The fact that the Hartree approximation 
does not satisfy Goldstone's theorem means that the resulting particles must be massive, 
explaining why somewhat fewer particles are produced. 

Finally, we show the special case in which there is no initial biasing in the field, r](to = 
0) = 0, rj(t o = 0) = 0, and H(t ) = 0.1 in figures 3a-d. With such an initial condition, the 
Hartree approximation and the large N limit are equivalent. The zero mode remains zero 
for all time, so that the quantity gH,(t) (fig. 3a) satisfies the sum rule ( |4.9|) by reaching the 
value one without decaying for late times. Notice that many more particles are produced in 
this case (fig 3b); the growth of the particle number for late times is due to the expansion 
of the universe. The particle distribution (fig. 3c) is similar to that of the slow roll case in 
figure 1. The equation of state (fig. 3d) is likewise similar. 

In each of these cases of slow roll dynamics, increasing the Hubble constant has the 
effect of slowing the growth of both rj and gH. The equation of state will be that of a de 
Sitter universe for a longer period before moving to a matter dominated equation of state. 
Otherwise, the dynamics is much the same as in figs. 1-3. 

Case 2: m 2 < 0, i](t ) ^> 1. We now examine the case of preheating occurring in a 
chaotic inflationary scenario with a symmetry broken potential. In chaotic inflation, the 
zero mode begins with a value rj(t) 3> 1. During the de Sitter phase, H ^> 1, and the 
field initially evolves classically, dominated by the first order derivative term appearing in 
the zero mode equation (see ( 2.11|) and ([2.21 )). Eventually, the zero mode rolls down the 



potential, ending the de Sitter phase and beginning the preheating phase. We consider the 
field dynamics in the FRW universe where preheating occurs after the end of inflation. We 
thus take the initial temperature to be zero, T = 0. 

Figure 4 shows our results for the quantities, r](t), gT,(t), gN(t), gN k (t), and p(t)/e(t) 
for the evolution in the large N approximation within a radiation dominated gravitational 
background with H(to) = 0.1. The initial condition on the zero mode is chosen to have the 
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representative value 77(^0 = 0) = 4 with r](to = 0) = 0. Initial values of the zero mode much 
smaller than this will not produce significant growth of quantum fluctuations; initial values 
larger than this produces qualitatively similar results, although the resulting number of 
particles will be greater and the time it takes for the zero mode to settle into its asymptotic 
state will be longer. 

We see from figure 4a that the zero mode oscillates rapidly, while the amplitude of 
the oscillation decreases due to the expansion of the universe. This oscillation induces 
particle production through the process of parametric amplification (fig. 4c) and causes the 
fluctuation gll to grow (fig. 4b). Eventually, the zero mode loses enough energy that it 
is restricted to one of the two minima of the tree level effective potential. The subsequent 
evolution closely follows that of Case 1 above with gll decaying in time as l/a 2 (t) ~ I ft 
with 7] given by the sum rule (|4.9|) . The spectrum (fig. 4d) indicates a single unstable band 
of particle production dominated by the modes k — 1/2 to about k = 3 for late times. The 
structure within this band becomes more complex with time and shifts somewhat toward 
lower momentum modes. Such a shift has already been observed in Minkowski spacetimes [|7| . 
Figure 4e shows the equation of state which we see to be somewhere between the relations 
for matter and radiation for times out as far as t = 400, but slowly moving to a matter 
equation of state. Since matter redshifts as l/a 3 (t) while radiation redshifts as l/a 4 (£), the 
equation of state should eventually become matter dominated. Given the equation of state 
indicated by fig. 4e, we estimate that this occurs for times of order t = 10 4 . The reason 
the equation of state in this case differs from that of cold matter as was seen in figs. 1-3 is 
that the particle distribution produced by parametric amplification is concentrated at higher 
momenta, k ~ 1. 

Figure 5 shows the corresponding case with a matter dominated background. The results 
are qualitatively very similar to those described for figure 4 above. Due to the faster expan- 
sion, the zero mode (fig. 5a) finds one of the two wells more quickly and slightly less particles 
are produced. For late times, the fluctuation g£ (fig. 5b) decays as l/a 2 (t) oc 1/t 4 / 3 . Again 
we see an equation of state (figs. 5e) which evolves from a state between that of pure 
radiation or matter toward one of cold matter. 

The Hartree case is depicted in figure 6 for a matter dominated universe, with the initial 
condition on the zero mode 77(^0 = 0) = 4y / 3- Again, the evolution begins in much the same 
manner as in the large iV approximation with oscillation of the zero mode (fig. 6a), which 
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eventually settles into one of the two minima of the effective potential. Whereas in the 
large N approximation, the zero mode approaches the minimum asymptotically [as given by 
( |4.9|) and our late time analysis below], in the Hartree approximation we see that the zero 
mode finds the minimum quickly and proceeds to oscillate about that value. The two point 
correlator (fig. 6b) quickly evolves toward zero without growing large. Particle production 
in the Hartree approximation (figs. 6c-d) is again seen to be inefficient compared to that of 
the large N case above. Fig. 6e again shows that the equation of state is matter dominated 
for all but the earliest times. 

A larger Hubble constant prevents significant particle production unless the initial ampli- 
tude of the zero mode is likewise increased such that the relation 77 (to) 3> H(to) is satisfied. 
For very large amplitude 77 (to) ^ 1 5 to the extent that the mass term can be neglected 
and while the quantum fluctuation term has not grown to be large, the equations of motion 
( |2.11| ), ( |2.12| ), and (2.21) are scale invariant with the scaling rj —>■ firj, H —>■ fiH, t —>■ t/fi, 
and k — > /ik, where /x is an arbitrary scale. 

For completeness, we show the case of the evolution with initial values of the Hubble 
constant given by H(t ) = 5 and H(t ) = 2 respectively in figures 7-8 using radiation 
dominated expansion. Here, we have used the large iV approximation and have made an 
appropriate increase in the initial value of the zero mode such that the fluctuations (figs. 
7b and 8b) grow significantly [we have chosen rj(t ) = 40 in fig. 7 and 7/(t ) = 16 in fig. 8]. 
While the dynamics looks much like that of figs. 4-6 above, we point out that the particle 
distribution (figs. 7d and 8d) is extended to higher values with the result being that the 
equation of state (figs. 7e and 8e) is weighted more toward that of radiation. 

Case 3: 771 2 > 0, 7?(to) 3> 1. The final case we examine is that of a simple chaotic scenario 
with a positive mass term in the Lagrangian. Again, preheating can begin only after the 
inflationary phase of exponential expansion; this allows us to take a zero temperature initial 
state for the FRW stage. 

Figure 9 shows this situation in the large iV approximation for a matter dominated 
cosmology. The zero mode, rj(t), oscillates in time while decaying in amplitude from its 
initial value of 7/(to = 0) = 5, rjito = 0) = (fig. 9a), while the quantum fluctuation, <?E, 
grows rapidly for early times due to parametric resonance (figs. 9b). We choose here an 
initial condition on the zero mode which differs from that of figs 4-5 above since there is no 
significant growth of quantum fluctuations for smaller initial values. From figure 9d, we see 
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that there exists a single unstable band at values of roughly k = 1 to k = 3, although careful 
examination reveals that the unstable band extends all the way to k = 0. The equation of 
state is depicted by the quantity p(t)/e(t) in figure 9e. As expected in this massive theory, 
the equation of state is matter dominated. 

The final case is the Hartree approximation, shown in figure 10. Here, parametric am- 
plification is entirely inefficient when expansion of the universe is included and we require 
an initial condition on the zero mode of 77 (to = 0) = 12\/3 to provide even meager growth 
of quantum fluctuations. We have used a matter dominated gravitational background with 
H(t ) =0.1. We see that while the zero mode oscillates (fig. 10a), there is little growth 
in quantum fluctuations (fig. 10b) and few particles produced (fig. 10c). Examining the 
particle distribution (fig. lOd), it is found that the bulk of these particles is produced within 
a single resonance band extending from fc~15toA;~16. This resonance develops at early 
time during the large amplitude oscillation of the zero mode. These results are explained 
by a simple resonance band analysis described below. 

At first glace, it is not entirely clear why there are so many more particles produced in 
the large N case of figure 9 than in the Hartree case of figure 10. Since in the present case 
the Hubble time is long compared to the oscillation timescale of the zero mode, H <C 1, 
we would expect a forbidden band for early times at the location given approximately by 
the Minkowski results provided in Ref. 0. In fact, we find this to be the case. However, 
we know from previous studies that in Minkowski space a similar number of particles is 
produced in both the Hartree and large N case HH. 

The solution to this problem is inherent in the band structure of the two cases when 
combined with an understanding of the dynamics in an expanding spacetime. First, we note 
that, for early times when g£ <C 1, the zero mode is well fit by the function r)(t) = t]of{t) / a(t) 
where f(t) is an oscillatory function taking on values from —1 to 1. This is clearly seen from 
the envelope function rjo/a(t) shown in fig. 10a (recall that gH <C 1 during the entire 
evolution in this case). Second, the momentum that appears in the equations for the modes 
( p.!2| ) is the physical momentum k/a(t). We therefore write the approximate expressions 
for the locations of the forbidden bands in FRW by using the Minkowski results of || with 
the substitutions t/q — > 77?Q/a 2 (t) (where the factor of 7 accounts for the difference in the 
definition of the non-linear coupling between this study and ||) and q 2 — > k 2 /a 2 (t). 

Making these substitutions, we find for the location in comoving momentum k of the 
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forbidden band in the large N (fig. 9) and Hartree (fig. 10) cases: 

2 

< k 2 < 5°, (largeN) (4.11) 



^ + 3a 2 (t) < k 2 <a 2 (t) 



^ + ^TT + 4 + 1 I • {Hartree) (4.12) 



2 w " - y 1 \ \ 3a 4 (t) a 2 (t) 

The important feature to notice is that while the location of the unstable band (to a first 
approximation) in the case of the continuous O(N) theory is the same as in Minkowski and 
does not change in time, the location of the band is time dependent in the discrete theory 
described by the non-perturbative Hartree approximation. 
While r] /a(t) ^> 1, the Hartree relation reduces to 

"A < k 2 < 4=- (4.13) 
2 v3 

This is the same as the Minkowski result for large amplitude, and one finds that this expres- 
sion accurately predicts the location of the resonance band of fig. lOd. However, with time 
the band shifts its location toward higher values of comoving momentum as given by fl4.12Q , 
cutting off particle production in that initial band. There is continuing particle production 
for higher modes, but since the Floquet index is decreased due to the reduced amplitude 
of the zero mode, since there is no enhancement of production of particles in these modes 
(as these modes begin with at most of order 1 particles), and because the band continues 
to shift to higher momenta while becoming smaller in width, this particle production never 
becomes significant. 

As in the symmetry broken case of figs. 4-6, the equations of motion for large amplitude 
and relatively early times are approximately scale invariant. In figures 11-12 we show the 
case of the large N evolution in a radiation dominated universe with initial Hubble constants 
of H(to) = 5 and H(to) = 2 respectively with appropriately scaled initial values of the 
zero mode of 77(^0) = 40 and 7?(£o) = 16. Again, the qualitative dynamics remains largely 
unchanged from the case of a smaller Hubble constant. 

C. Late Time Behavior 

We see clearly from the numerical evolution that in the case of a symmetry broken 
potential, the late time large N solutions obey the sum rule 



20 



0. 



(4.14) 



This sum rule is a consequence of the late time Ward identities which enforce Goldstone's 
Theorem. Because of this sum rule, we can write down the analytical expressions for the 
late time behavior of the fluctuations and the zero mode. Using Q4.14| ), the mode equation 
( [2.12j ) becomes 

k 2 



d 2 a(t) d 
d^ + d~e)Jt + a 2 (t) 



U k (t) = 0. 



(4.15) 



This equation can be solved exactly if we assume a power law dependence for the scale factor 
a(t) = (t/t ) n . The solution is 



U k (t) = c k t^-^' 2 



2-2n 



n — 1 



+ d k t^-^' 2 Yu 



'kt^t 1 -^ 
n — 1 



(4.16) 



where J u and Y v are Bessel and Neumann functions respectively, and the constants c k and 
d k carry dependence on the initial conditions and on the dynamics up to the point at which 
the sum rule is satisfied. 

These functions have several important properties. In particular, in radiation or matter 
dominated universes, n < 1, and for values of wavenumber satisfying k t~^~ n '/tQ, the 
mode functions decay in time as l/a(t) ~ t~ n . Since the sum rule applies for late times, 
t — to ^> 1 in dimensionless units, we see that all values of k except a very small band about 
k = redshift as l/a(t). The k = mode, however, remains constant in time, explaining 
the support evidenced in the numerical results for values of small k (see figs. 1,3). These 
results mean that the quantum fluctuation has a late time dependence of (tp 2 (t)) r ~ l/a 2 (t). 
The late time dependence of the zero mode is given by this expression combined with the 
sum rule ( f4.14| ). These results are accurately reproduced by our numerical analysis. Note 
that qualitatively this late time dependence is independent of the choice of initial conditions 
for the zero mode, except that there is no growth of modes near k = in the case in which 
particles are produced via parametric amplification (figs. 4,5). 

For the radiation n = ~ and matter dominated n = | universes, eq.( [4.16| ) reduces to 
elementary functions: 



a(t) U k (t) = c k e 2ikt o 



2iktl /2 t^ 



+ d k e 



(RD) , 



a(t) U k (t) = c k e 3lkt ° 3tl/3 



1 + 



ZktTt 1 '^ 



+ d k e 6 o 



3Ho /3 t 1/3 . 



(MD) . (4.17) 
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It is also of interest to examine the n > 1 case. Here, the modes of interest satisfy the 
condition k <C t n_1 /tg for late times. These modes are constant in time and one sees that 
the modes are frozen. In the case of a de Sitter universe, we can formally take the limit 
n — > oo and we see that all modes become frozen at late times. This case was studied in 
detail in ref. flTSJ. 

V. CONCLUSIONS 

We have shown that there can be significant particle production through quantum fluc- 
tuations after inflation. However, this production is somewhat sensitive to the expansion 
of the universe. From our analysis of the equation of state, we see that the late time dy- 
namics is given by a matter dominated cosmology. We have also shown that the quantum 
fluctuations of the inflaton decay for late times as l/a 2 (t), while in the case of a symmetry 
broken inflationary model, the inflaton field moves to the minimum of its tree level poten- 
tial. The exception to this behavior is the case when the inflaton begins exactly at the 
unstable extremum of its potential for which the fluctuations grow out to the minimum of 
the potential and do not decay. Initial production of particles due to parametric amplifi- 
cation is significantly greater in chaotic scenarios with symmetry broken potentials than in 
the corresponding theories with positive mass terms in the Lagrangian, given similar initial 
conditions on the zero mode of the inflaton. 

Since there are a number of articles in the literature treating the problem of preheating, 
it is useful to review the unique features of the present work. First, we have treated the 
problem dynamically, without using the effective potential (an equilibrium construct) to 
determine the evolution. Second, we have provided consistent non-perturbative calculations 
of the evolution to bring out some of the most relevant aspects of the late time behavior. In 
particular, we found that the quantum backreaction naturally inhibits catastrophic growth of 
fluctuations and provides a smooth transition to the late time regime in which the quantum 
fluctuations decay as the zero mode approaches its asymptotic state. Third, the dynamics 
studied obeys the constraint of covariant conservation of the energy momentum tensor. 

The next stage of this analysis is to allow dynamical evolution of the scale factor, to 
follow the evolution of the inflaton through the de Sitter stage and into the stage of particle 
production. As the expansion rates for which there is significant particle production is some- 
what restrictive, it is yet to be seen whether in a theory with fully dynamical gravitational 
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expansion such particle production will be a significant factor. This analysis is currently 
underway [34 . 
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APPENDIX A: CONFORMAL TIME ANALYSIS AND INITIAL CONDITIONS 

The issue of renormalization and initial conditions is best understood in conformal time 
which is a natural framework for adiabatic renormalization and regularization. 
Quantization in conformal time proceeds by writing the metric element as 

ds 2 = C 2 (r)(dT 2 - dx 2 ). (Al) 

Under a conformal rescaling of the field 

^(x,t)=x(x,T)/C(r), (A2) 

the action for a scalar field (with the obvious generalization to N components) becomes, 
after an integration by parts and dropping a surface term 

S = J d'xdr |i( X ') 2 - \{V X f - V(x)} , (A3) 

with 

V( X ) = C\r)V^/C{r)) - C 2 (r)| X 2 , (A4) 

where 1Z = QC"(t)/C 3 (t) is the Ricci scalar, and primes stand for derivatives with respect 
to conformal time r. 

The conformal time Hamiltonian operator, which is the generator of translations in r, is 
given by 

H T = J d 3 x [\u 2 x + l -{V X f + V(x)} , (A5) 
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with II X being the canonical momentum conjugate to x-, n x = x' ■ Separating the zero mode 



of the field x 



x(x,t) = xo(t) + x(x,t), 



(A6) 



and performing the large N or Hartree factorization on the fluctuations we find that the 
Hamiltonian becomes linear plus quadratic in the fluctuations, and similar to a Minkowski 
space-time Hamiltonian with a r dependent mass term given by 



M 2 (t) 



C 2 (t) 



™ 2 + tt-\)n + \xl{r) + \(x 2 ) 



(A7) 



We can now follow the steps and use the results of reference [T7| for the conformal time 
evolution of the density matrix by setting a(t) = 1 in the proper equations of that reference 
and replacing the frequencies by 



ul{r) = k 2 + M\r), 



(A8) 



and the expectation value in ( |AT| ) is obtained in this r evolved density matrix. The time 
evolution of the kernels in the density matrix (see |17|]) is determined by the mode functions 
that obey 



dr 1 



fk{r) = 0. 



(A9) 



The Wronskian of these mode functions 



/*) = f'Jl - fkfk 



(A10) 



is a constant. It is natural to impose initial conditions such that at the initial r the density 
matrix describes a situation of local thermodynamic equilibrium and therefore commutes 
with the conformal time Hamiltonian at the initial time. This implies that the initial con- 
ditions of the mode functions fkij) should be chosen to be (see [[17]]) 

1 



fk(r ) 



fk( T o) = -i\]^k{r )fk{T ). 



(All) 



These initial conditions correspond to the choice of mode functions which coincide with the 



first order adiabatic modes and those of the Bunch-Davies vacuum for large momentum |29 



To see this clearly, consider the WKB solutions of the mode equation (|A9|) of the form 
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with the function R k {r) obeying the Riccati equation 

R' k + R 2 k + k 2 + M 2 (t) =0. 



(A12) 



(A13) 



This equation possesses the solution 

#fc( r ) = -tk + R ,k{ T ) ~ 1 — ; 1 — 1 — — 1 1- ' 



k k 2 k 3 

and its complex conjugate. We find for the coefficients: 

1 



A: 4 



Ro,k = ; R 1)k = -M\t) ; R 2 , k 



2 ^ 



y ( R 3,k + 2-Ri,fc-R 2 ,, 



(A14) 



(A15) 



The solutions /a;(t) obeying the boundary conditions ([All]) are obtained as linear combina- 
tions of this WKB solution and its complex conjugate 

1 



fk(r) 



2yW(To) 



[(l + 7 )L> fc ( T ) + (l- 7 )L>*( T )] 



(A16) 



where the coefficient 7 is obtained from the initial conditions. It is straightforward to find 
that the real and imaginary parts are given by 



7r = 1 + 0(1A 4 ) ; ll = 0{l/k 3 ). 



(A17) 



Therefore the large-/c mode functions satisfy the adiabatic vacuum initial conditions |29 
This, in fact, is the rationale for the choice of the initial conditions (All 



Following the analysis presented in [17] we find, in conformal time that 

d 3 k 



\fk(r)\" 



(A18) 



2(2vr) 3 

The Heisenberg field operators x(x, r) and their canonical momenta U x (x,t) can be ex- 
panded as: 

d 3 k 



X(x, t) 



n x (£, r) 



x/2(2tt) 3 / 2 
d 3 k 
v / 2(2tt) 3 / 2 



ik-x 



a k f k (r) + aLJt (r 



ik-x 



(A19) 
(A20) 
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with the time independent creation and annihilation operators and a\ obeying canonical 
commutation relations. Since the fluctuation fields in comoving and conformal time are 
related by a conformal rescaling 

X(x, t) 



ip(x, t) 



C(r) 



(A21) 



it is straightforward to see that the mode functions in comoving time are related to those in 
conformal time simply as 

fk(r) 



U k (t) 



C(r) 



(A22) 



Therefore the initial conditions ( |A1 1|) on the conformal time mode functions imply the initial 
conditions for the mode functions in comoving time are given by 

1 



£4 (to) 



U k (h 



-iuJk(T ) - H(t )] U k (to), 



(A23) 



where we have chosen the normalization of the scale factor such that a(t ) = C(t ) = 1. 

For renormalization purposes we need the large- k behavior of \Uk{t)\ 2 , \ Uk{t)\ 2 , which are 
determined by the large-fc behavior of the conformal time mode functions and its derivative. 
These are given by 



l/*(r)| s 



RxAj) 1 
k 2 k 4 



'R!{ 



k\ T ) , 3 2 / \ | 

+ n R lA T ) + 



4 



(A24) 



\f'k(r) 



k 



1 + + F T~ + 2^ fc(r) + 



(A25) 



We note that the large k behavior of the mode functions to the order needed to renor- 
malize the quadratic and logarithmic divergences is insensitive to the initial conditions. This 
situation must be contrasted with the case in which the initial conditions in comoving time 



are imposed as described in Thus the merit in considering the initial conditions in 

conformal time described in this article. 

The correspondence with the comoving time mode functions is given by: 

\h{r)\ 2 



\Uu(t)\ 2 

W)\ 2 



C 2 (r) 



1 



C 2 



T 



C 2 (r) 



H 2 



H d ' 



IA(r)| s 



(A26) 



These are asymptotic forms used in the renormalization program in section III. 
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There is an important physical consequence of this choice of initial conditions, which is 
revealed by analyzing the evolution of the density matrix. 

In the large iV or Hartree (also to one-loop) approximation, the density matrix is Gaus- 
sian, and defined by a normalization factor, a complex covariance that determines the diag- 
onal matrix elements and a real covariance that determines the mixing in the Schrodinger 
representation as discussed in reference [|17|] (and references therein). 

In conformal time quantization and in the Schrodinger representation in which the field 
X is diagonal the conformal time evolution of the density matrix is via the conformal time 
Hamiltonian (|A5|). The evolution equations for the covariances is obtained from those given 
in reference [DJ by setting a(t) = 1 and using the frequencies o;|(r) = k 2 + A4 2 (r). In 
particular, by setting the covariance of the diagonal elements (given by equation (2.20) in 
jl7ll ; see also equation (2.44) of |jT7|), 

Mr) = _i 4?TT' (A27) 
Ik v) 

we find that with the initial conditions (|A11|) , the conformal time density matrix is that 
of local equilibrium at tq in the sense that it commutes with the conformal time Hamilto- 
nian. However, it is straightforward to see, that the comoving time density matrix does not 
commute with the comoving time Hamiltonian at the initial time to- 

An important corollary of this analysis and comparison with other initial conditions used 
in comoving time is that assuming initial conditions of local equilibrium in comoving time 
leads to divergences that depend on the initial condition as discussed at length in [17]. This 
dependence of the renormalization counterterms on the initial condition was also realized by 
Leutwyler and Mallik |2l| within the context of the CTP formulation. Imposing the initial 
conditions corresponding to local thermal equilibrium in conformal time, we see that: i) 
the renormalization counterterms do not depend on the initial conditions and ii) the mode 
functions are identified with those corresponding to the adiabatic vacuum for large momenta. 

Thus this analysis justifies the use of the initial conditions on the comoving mode func- 
tions ( |A23|) . Furthermore, being that the comoving time density matrix does not describe at 
any time a condition of local thermodynamic equilibrium, the "temperature" that enters in 
the mixing covariance in the density matrix is understood as a parameter describing a mixed 
state with the notion of temperature at local thermodynamic equilibrium in conformal time 
quantization. 
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For our main analysis we choose this "temperature" to be zero so that the resulting 
density matrix describes a pure state, which for the large momentum modes coincides with 
the conformal adiabatic vacuum. 

Particle Number: 

We write the Fourier components of the field \ and its canonical momentum Il x given 
by flftiq ) - (|X20D as: 

1 

V2 



Xk{T) = 

n x ,fcO-) 



a k fk{r) +al k fl(r) 



1 r 



(A28) 
(A29) 



These (conformal time) Heisenberg operators can be written equivalently in terms of the r 
dependent creation and annihilation operators 

, , 1 



Xk{T) 



a fc (r)e-^ (ro)r + ~al(T)e luJk(To)T 



n 



X,k\ T ) 



a k {r)e 



-iuj k (T )T _ ~f 



ai(T e 



(A30) 



(A31) 



The operators a k (r) ; a k (r) are related by a Bogoliubov transformation. The number of 
particles referred to the initial Fock vacuum of the modes f k , is given by 



NJr) 



'a k (r)a k (r)) 



fk(r) 



fk(0) 



1 + 



Wfc(ro) 



fL(r) 



fk(r) 



(A32) 



or alternatively, in terms of the comoving mode functions U k (t) = f k (T)/C(r) we find 



N k (t) 



a 2 (t) 



U k {t) 



£4(0) 



l + 



^(0) 



U k {t) + HU k {t) 



U k (t) 



1 
2 



(A33) 



Using the large /c-expansion of the conformal mode functions given by eqns. (|A24j ) -( A25 ) 
we find the large-fc behavior of the particle number to be N k fc =°° 0(1/ k 4 ), and the total 
number of particles (with reference to the initial state at r ) is therefore finite. 
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FIGURES 
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FIG. 1. Figure 1: Symmetry broken, slow roll, large N, matter dominated evolution of (a) the 
zero mode 77(f) vs. t, (b) the quantum fluctuation operator gT,(t) vs. t, (c) the number of particles 
gN(t) vs. t, (d) the particle distribution gNk(t) vs. k at t = 149.1 (dashed line) and t = 398.2 
(solid line), and (e) the ratio of the pressure and energy density p(t)/e(t) vs. t for the parameter 
values m 2 = -1, r/(t ) = 10" 7 , r)(t ) = 0, g = 1(T 12 , H(t ) = 0.1. 
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Figure 2a 
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FIG. 2. Figure 2: Symmetry broken, slow roll, Hartree, matter dominated evolution of (a) the 
zero mode r](t) vs. t, (b) the quantum fluctuation operator gT,(t) vs. t, (c) the number of particles 
gN(t) vs. t, (d) the particle distribution gNk{t) vs. k at t = 150.7 (dashed line) and t = 396.1 
(solid line), and (e) the ratio of the pressure and energy density p(t)/e(t) vs. t for the parameter 
values m 2 = -1, r](t ) = 3 1 / 2 • 10~ 7 , r}{t ) = 0, g = 1(T 12 , H(t ) = 0.1. 
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Figure 3a 



Figure 3b 
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FIG. 3. Figure 3: Symmetry broken, no roll, matter dominated evolution of (a) the quantum 
fluctuation operator gT>(t) vs. t, (b) the number of particles gN(t) vs. t, (c) the particle distribution 
gNk(t) vs. k at t = 150.1 (dashed line) and t = 397.1 (solid line), and (d) the ratio of the pressure 
and energy density p(t)/e(t) vs. t for the parameter values m 2 = —1, r/(io) = 0, r)(to) = 0, 
g= 10- 12 , H(t ) = 0.1. 
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FIG. 4. Figure 4: Symmetry broken, chaotic, large N, radiation dominated evolution of (a) the 
zero mode rj(t) vs. t, (b) the quantum fluctuation operator gT,(t) vs. t, (c) the number of particles 
gN(t) vs. t, (d) the particle distribution gNk(t) vs. k at t = 76.4 (dashed line) and t = 392.8 (solid 
line), and (e) the ratio of the pressure and energy density p(t)/e(t) vs. t for the parameter values 
m 2 = -1, 7/(t ) = 4, j)(t ) = 0, g = 1(T 12 , fl"(t ) = 0.1. 
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Figure 5a 
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FIG. 5. Figure 5: Symmetry broken, chaotic, large N, matter dominated evolution of (a) the 
zero mode rj(t) vs. t, (b) the quantum fluctuation operator gT,(t) vs. t, (c) the number of particles 
gN(t) vs. t, (d) the particle distribution gNk{t) vs. k at t = 50.8 (dashed line) and t = 399.4 (solid 
line), and (e) the ratio of the pressure and energy density p(t)/e(t) vs. t for the parameter values 



m 



-1, r,(t ) = 4, j)(t ) = o, 5 = irr 12 , fl-(to) = 0.1. 
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Figure 6a Figure 6b 
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FIG. 6. Figure 6: Symmetry broken, chaotic, Hartree, matter dominated evolution of (a) the 
zero mode rj(t) vs. t, (b) the quantum fluctuation operator gT,(t) vs. t, (c) the number of particles 
gN(t) vs. t, (d) the particle distribution gNk(t) vs. k at t = 151.3 (dashed line) and t = 397.0 
(solid line), and (e) the ratio of the pressure and energy density p(t)/e(t) vs. t for the parameter 
values m 2 = -1, r](t ) = 4 • 3 1 / 2 , rj{t ) = 0, g = 1(T 12 , H(t ) = 0.1. 
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Figure 7a 
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FIG. 7. Figure 7: Symmetry broken, chaotic, large N, radiation dominated evolution of (a) the 
zero mode rj(t) vs. t, (b) the quantum fluctuation operator <?£(t) vs. t, (c) the number of particles 
gN(t) vs. t, (d) the particle distribution gNk(t) vs. k at t = 118.9 (dashed line) and t = 394.7 
(solid line), and (e) the ratio of the pressure and energy density p(t)/e(t) vs. t for the parameter 
values m 2 = -1, r}(t ) = 40, r](t ) = 0, g = 1(T 12 , H(t ) = 5.0. 
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FIG. 8. Figure 8: Symmetry broken, chaotic, large N, radiation dominated evolution of (a) the 
zero mode rj(t) vs. t, (b) the quantum fluctuation operator gT,(t) vs. t, (c) the number of particles 
gN(t) vs. t, (d) the particle distribution gNk{t) vs. k at t = 55.1 (dashed line) and t = 194.2 (solid 
line), and (e) the ratio of the pressure and energy density p(t)/e(t) vs. t for the parameter values 



m 



•1, r,(t ) = 16, r,(to) =0,g= 10" 12 , H(t ) = 2.0 
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Figure 9a 
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FIG. 9. Figure 9: Symmetry unbroken, chaotic, large N, matter dominated evolution of (a) the 
zero mode rj(t) vs. t, (b) the quantum fluctuation operator <?£(t) vs. t, (c) the number of particles 
gN(t) vs. t, (d) the particle distribution gNk{t) vs. k at t = 77. 4 (dashed line) and t = 399.7 (solid 
line), and (e) the ratio of the pressure and energy density p(t)/e(t) vs. t for the parameter values 
m 2 = +1, rj(t ) = 5, j)(t ) = 0, g = KT 12 , H(t ) = 0.1. 
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Figure 10a 
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Figure 10b 
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FIG. 10. Figure 10: Symmetry unbroken, chaotic, Hartree, matter dominated evolution of (a) 
the zero mode r)(t) vs. t, (b) the quantum fluctuation operator <?£(i) vs. t, (c) the number of 
particles gN(t) vs. t, (d) the particle distribution gNk(t) vs. k at t = 50.5 (dashed line) and 
t = 391.2 (solid line), and (e) the ratio of the pressure and energy density p(t)/e(t) vs. t for the 
parameter values m 2 = +1, r](t ) = 12 • 3 1 / 2 , fj(t ) = 0, g = 10~ 12 , H(t ) = 0.1. 
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FIG. 11. Figure 11: Symmetry unbroken, chaotic, large TV, radiation dominated evolution of 
(a) the zero mode r)(t) vs. t, (b) the quantum fluctuation operator <?£(i) vs. t, (c) the number 
of particles gN(t) vs. t, (d) the particle distribution gNk{t) vs. k at t = 117.3 (dashed line) and 
t = 393.6 (solid line), and (e) the ratio of the pressure and energy density p(t)/e(t) vs. t for the 
parameter values m 2 = +1, r](t ) = 40, f](t ) = 0, g = 1CT 12 , H(t ) = 5.0. 
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Figure 12a 
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FIG. 12. Figure 12: Symmetry unbroken, chaotic, large TV, radiation dominated evolution of 
(a) the zero mode r\{t) vs. t, (b) the quantum fluctuation operator <?£(i) vs. t, (c) the number 
of particles gN(t) vs. t, (d) the particle distribution gNk{t) vs. k at t = 102.1 (dashed line) and 
t = 251.6 (solid line), and (e) the ratio of the pressure and energy density p(t)/e(t) vs. t for the 
parameter values m 2 = +1, r](t ) = 16, f](t ) = 0, g = 1CT 12 , H(t ) = 2.0. 



43 



